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It is shown for an electron-nucleus mixture that the electron and nuclear pressures are 
defined clearly and simply by the virial theorem; the total pressure of this system is a sum of 
these two pressures. The electron pressure is different from the conventional electron pressure 
being expressed as the sum of two times of kinetic energy and the potential energy in that 
the nuclear virial term is subtracted; this fact is exemplified by several kinds of definitions for 
the electron pressure enumerated in this work. The conventional definition of the electron 
pressure in terms of the nuclear virial term is shown inappropriate. Similar remarks are 
made about the definition of the stress tensor in this mixture. It is also demonstrated that 
both of the electron and nuclear pressures become zero at the same time for a metal in the 
vacuum, in contrast to the conventional viewpoint that the zero pressure is realized by a 
result of the cancellation between the electron and nuclear pressures, each of which is not 
zero. On the basis of these facts, a simple equation of states for liquid metals is derived, and 
examined numerically for liquid alkaline metals by use of the quantum hypernetted chain 
equation and the Ashcroft model potential. 

§1. Introduction 

In order to determine the properties of solids, liquid metals and plasmas as 
electron-nucleus mixtures, the nuclei (ions) are usually treated as classical particles, 
which are considered as giving an external potential for the electrons: as a result, 
we can apply the density-functional (DF) theory to this inhomogeneous electrons. 
Thus, the DF theory provides the inhomogeneous properties of the electrons under 
this external potential. Therefore, some problems occur to determine homogeneous 
quantities of an electron-nucleus mixture including the nuclear contribution from 
the results of the DF theory, where the nuclei appear only as providing an external 
potential. Such an example of problems occurs in the determination of the total 
pressure of electron-nucleus mixtures on the basis of the DF theory. The electron 
pressure is frequently discussed from the DF theory; nevertheless, the term, 'electron 
pressure' P e , is differently used in two ways in its meaning. Some investigators 1 -*' 2 ) 
take the electron pressure P e for the electrons confined in a volume V together 
with the fixed nuclei at the positions {R Q } to be 3P C V = 2T e + U e in terms of the 
exact kinetic and potential energies, T e and U e , respectively. On the other hand, 
many kinds of definitions for the electron pressure lead to an expression 3P e V = 
2T e + U e — R a • F Q with use of the force F a on a-nucleus, as will be shown in this 
work. This difference is not only a problem of the definition, since the total pressure 
of an electron-nucleus mixture as a homogeneous system is to be determined from 
the electron pressure in conjunction with the nuclear contribution: the definition 
of electron pressure should be consistent with this purpose, as will be done in this 
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work. Therefore, it is important to see the meaning of the electron pressure in the 
calculation of the total pressure taking account of the nuclear contribution. 

Following Slater, 3 ) many investigators 1 )' 2 )' 4 ) have identified the electron pres- 
sure with the total force on the nuclei; 3P e V = J2 a ' ^a- As a consequence of 
this definition which leads to the relation 3P C V = 2T e + U e , the 'electron pressure' 
becomes different from others as mentioned above. Similarly, the same problem was 
seen in the definition of the stress tensors 4 ) in terms of the total force on the nu- 
clei. This definition of P e is based on the assumption: an electron-nucleus mixture 
confined in a volume V by wall potentials, U^ e and U^ nc , for the electrons and the 
nuclei, respectively, is thermodynamically identical with the mixture with only the 
nuclei confined in V by U™ c , since in the latter system the electrons seem to spill 
out outside the nuclei a distance only on an atomic scale even if U^ e = 0. In this 
work, we will examine the validity of this assumption with use of the virial theorem 
which leads to a definition of the electron pressure in a more general form. From 
this new definition of the electron pressure, the total pressure of an electron-nucleus 
mixture is shown to be represented as a sum of the electron and nuclear pressures. 

Also, we can give an important remark about the fact that the pressure in 
solids or liquid metals becomes zero in the vacuum; it is a standard viewpoint 5 )' 6 ) 
that the electron pressure is negative because of the electrons yielding a main part of 
cohesive energy, while the nuclear pressure is positive, and the total pressure becomes 
zero as a result of cancellation between two pressures in the vacuum. Contrary to 
this assertion, the total pressure becomes zero, since both the electron and nuclear 
pressures are zero at the same time in the vacuum, as we will see. In this work, we 
discuss on these problems from fundamental and simple virial relations. 

Based on the fact that the total pressure of an electron-nucleus mixture is a sum 
of the electron and nuclear pressures, we obtain a simpler pressure formula for a 
liquid metal, where the effective interaction among nuclei (ions) are approximately 
represented by a pair interaction; the nuclear pressure is expressed in the virial 
form using this pair interaction and the radial distribution function. On the other 
hand, the electron pressure in a simple metal may be described by the jellium model. 
However, the electron pressure determined from the jellium model becomes zero only 
at the electron density of r s ~ 4 and the jellium surface energy becomes negative 
for r s ~ 2, while the jellium bulk modulus is negative at r s ~ 6. These drawbacks 
of the jellium model are rectified by the stabilized jellium (SJ) model 7 ) (or the ideal 
metal model 8 )' 9 )), where an infinitesimally thin dipole layer is added to the surface 
of the uniform background in the jellium; this dipole layer produces a uniform field 
in the jellium to make it stable at any metallic densities. On use of this SJ model 
to obtain an electron pressure expression, we derived a pressure formula for a simple 
liquid metal in the present work. 

In the next section, we show that the electron and nuclear pressures are clearly 
defined on the basis of the simple and fundamental virial theorem, and construct 
three systems by changing the wall potentials confining the electrons and the nuclei 
in a finite volume. Based on these systems, we make several statements about 
the pressures for an electron-nucleus mixture. Similar remarks are made about the 
definition of the stress tensors for this mixture in Sj3j In £JH we enumerate several 
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kinds of the definitions of the electron pressures proposed up to date. Equation of 
states (EOS) for liquid metals is set up in £15.11 a numerical examination of this EOS 
is performed in £15.21 The last section is devoted to summary and discussion. In 
Appendix [AJ the electron pressure formula is derived from the internal energy by 
performing a volume derivative in the framework of the DF theory. 

§2. Virial theorem for an electron-nucleus mixture 

The Hamiltonian for an electron-nucleus mixture is represented by 

^ = E^ + E^f + ^}'^})^^ + r N + f /, (2-i) 



where the coordinates and momenta of the electrons are denoted by {rj,pj} and 
those of the nuclei by {R Q ,P Q } with the Coulomb interactions between particles 



g-lV-J^-V Z * +iV , ZV . . (2.2) 
9 ^ If, — f,l ^ la. -a I 9^115 t> . I 



This Coulomb potential U of an electron-nucleus mixture satisfies the following op- 
erator relation 

U= -J^fi- ViL>-^R Q -V a J7. (2-3) 

Afterwards when the nuclei are regarded as classical particles, we use the symbols 
{Ro,,P a } instead of {R a ,P a }. In the standard approach to homogeneous liquid 
and solid metals, it is commonly assumed that the nuclei are treated as classical 
particles (or in the adiabatic approximation), and the electrons are considered as 
an inhomogeneous electron system under the external potential U ex caused by the 
nuclei fixed at their equilibrium positions. Here, it is important in treating an 
electron-nucleus mixture to recognize that there are four stages of approximations 
to the nuclear behaviour: (i) the nuclear mass is approximated as infinite, where the 
external potential U ex is determined uniquely from a nuclear configuration, since the 
nuclei do not move but remain to stay at their equilibrium positions; (ii) the nuclei 
are approximated as classical particles in the sense that its position and momentum 
are commutable; (iii) the nuclei are regarded as quantum particles, but can be treated 
in the adiabatic approximation; and finally (iv) the nuclei can not be treated in the 
adiabatic approximation. Therefore, we can apply the DF theory to the electrons 
in the electron-nucleus mixture without any approximation, if the nuclear mass is 
taken to be infinite. 

Generally, the thermodynamic pressure P of an inhomogeneous system caused 
by an external potential U ex (r) is clearly defined in terms of free energy F or internal 
energy E, 



dF 
dV 



IWex(r) 



dE 
dV 



(2-4) 

SNU ex (r) 
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under the condition that the volume parameter V of particles confined by a wall 
potential U w (r) remains to specify the state as a natural variable when the external 
potential U ex (r) is imposed, as is described in detail. 10 ^ In addition to the thermo- 
dynamic pressure, there is another 'kinetic' definition of the pressure using a wall 
potential U w on the basis of a standard assumption 11 )' 12 ) concerning the relation 
between the force exerted by the wall U w on the various particles and a hydrostatic 
pressure P: 

yVr* • Vit/w) = / Pr-dS = P<f r-dS = 3PV. (2-5) 
i JdV JdV 

Furthermore, the thermodynamic pressure of this system is shown to be identical 
with the 'kinetic' pressure in the form: 



dF 

3PV = -3V — 
oV 



= 2(f) - ■ V ^p) = • V ^w> , (2-6) 

TNU ex (v) i i 



where the potential operator U p is interpreted as involving the external potential 
C/e X , and T is the kinetic operator. 12 )~ 15 ) Therefore, it is necessary from (|2-5[) that 
the pressure takes a constant value at every point of the surface dV for the volume 
parameter V to be a natural variable. A simple example of this inhomogeneous 
system is an atom confined in a spherical wall with a nucleus fixed at its center. 

[ i-a ] As the first step to treat an electron-nucleus mixture on the basis of (|2-4j) 
which afford to use the DF theory, let us consider a mixture of electrons and nuclei 
with the infinite nuclear mass, confined in a finite volume V by wall potentials, U^ e 
and U^ nc , respectively. Here, we assume this mixture constitutes a crystalline state 
in equilibrium, which is used as a model for the electron-band calculations. Since 
the nuclei remain at their fixed equilibrium positions, the nuclei make an external 
potential U ex (r) for the electrons, and the volume V is a natural variable of this 
system. As a result, the pressure of this inhomogeneous electrons is determined by 
(!2-4p in terms of the electron free energy J- e , which is explicitly defined later in fl4-2f) . 
It is important to notice that the derivatives must be performed under the external 
potential being fixed. Therefore, the nuclei producing the external potential C/ ex for 
the electrons do not contribute to the pressure; the pressure is the electron pressure 
in the similar meaning to that of a noninteracting electron gas under the external 
potential. Since the electron free energy J- e can be described exactly in terms of the 
DF theory, 10 - 1 the electron pressure P e is shown (see Appendix |A|) to be 

= [2% + U e -J2^a- F a ]/W (2-7) 

TiVeC/ex(r) a 



dV 

[2(f e ) e - • V 4 f7) c ]/3V = [£<fi • X7 t U* c ) c ]/3V , (2 



where the bracket with a suffix e denotes the average over the electron states at a 
fixed nuclear configuration. Really, this can be also proven from (|2-6p and (|2-3p . 
Therefore, the electron pressure given above is a thermodynamic quantity, that is, 
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a physical quantity defined uniquely. Also, note that the electron pressure can be 
defined by a wall potential U^ e . However, it should be kept in mind that this pressure 
is not defined as a quantity of the homogeneous mixture. 

[ i-b ] In the second step to relate the above inhomogeneous system to a homo- 
geneous electron-nucleus mixture, we consider this system to be a binary mixture 
of electrons and nuclei. Then, the total pressure of this homogeneous mixture is 
written as 

3PV = 2(f c ) e - ■ V ^)c - • Vc ^> c 

i a 
i a 

= 3P C V + R«*-F Q = 2T C + U e = J2(ri ■ V^ e >e - Yl R -' F " ' ( 2 ' 9 ) 

a i a 

since the nuclei have infinite mass and do not move, but stay at their equilibrium 
positions. Note that the pressure of this homogeneous mixture is described in terms 
of the electron pressure P e defined by (12-7|) for the inhomogeneous system. Also, 
all of the nuclear forces F a are zero except for nuclei near the surface in the short 
range of the wall potential, since the nuclei are in their equilibrium positions. On 
the other hand, the wall forces satisfying = — F a near the surface describe the 
nuclear pressure. Thus, we can see that the total pressure is composed of the electron 
and nuclear pressures. In the derivation of (|2-9p . we have used the virial theorem for 
a homogeneous binary mixture confined in the volume V by a wall, which produces 
potentials, and U^, for A and B particles, respectively; 



dF 

3PV = -3V — 
ov 



TN 



i a 

^(r, • ViE#) + ]T<R a • V a C/£) . (2-10) 



[ ii ] Furthermore, when the nucleus has a finite mass, but behaves as a classical 
particle which makes the relations in the adiabatic approximation exact, the electron 
pressure of an homogeneous electron-nucleus mixture is given by the configuration 
average over {R Q } of the electron pressure in the inhomogeneous electrons under the 
external potential U cx (r\{R a }): 



3P C V=(-3V r&F ' 



dV 



) = 2{f e ) - Vft • V t U) , (2-11) 

T7V c (7 ox (r)/ N i 



as will be discussed later in f|6l [ iii ] Even when the nuclei behave as quantum par- 
ticles, Eq. (|2-lip is valid in the adiabatic approximation. 10 ) Therefore, this virial 
theorem (I2-I0j) suggests that we can define the electron pressure in a general form 
even when the adiabatic approximation becomes invalid, as is described in the fol- 
lowing. 
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[iv] From the Hamiltonian (|2-lj) . we can obtain the following relation: 

(Sy[fi • f>i,H]\& m ) = ih(V m \{fi/m - Vi ■ ViU}\* m ) =0, (2-12) 

due to the Hermitian property of the operator f j • p, of i-electron. Here, ^ m denotes 
an eigenfunction of the Hamiltonian H. The ensemble average of the electron kinetic 
energy is related to the interacting potential by use of (|2-12[) 

2<p?/2m) = (f< • ViU) . (2-13) 

In addition, we obtain the similar relation for a nucleus in the system 

2(P 2 /2M) = (R Q • V a U) . (2-14) 

Here, the ensemble average of an operator A is defined by 

(A) = Tr[exp(-/3tf)A]/Tr[exp(-/3#)] , (2-15) 

where the trace is taken over a complete set of states of the electron-nucleus mixture. 

Therefore, when the electrons and the nuclei are confined in a volume V by a 
wall with wall potentials, U^ e and U^ nc , respectively, Eqs. (|2-13|) and (12- 14h for this 
system are written in the following forms relating the electron and nuclear kinetic 
energies, T e and Tjy, respectively: 

2(f e ) = • Vfi) + • V^ le ) (2-16) 

i i 

2(Tn> = ^(R Q • Vjl) + ^(R Q • V a U™ c ) , (2-17) 

a a 

where the wall potentials should be taken into account in the ensemble average (|2- 15[) 
by adding them to the Hamiltonian 16 - 1 ' 17 ) or as the boundary condition 12 - 1-15 ) to solve 
the wave equation. At this point, it should be kept in mind that the interactions of 
the particles with the wall can be omitted from the new Hamiltonian if the boundary 
condition that W = for any particle on the wall is made. 12 )' 14 )' 15 ) As a result, with 
use of the Hamiltonian (|2- 1 j) and this boundary condition, the electron and nuclear 
pressures are expressed respectively in the forms: 

3P e V = • V,,C/; le ) = - J> • Ff) (2-18) 

i i 

= 2(f e ) - ■ Vi£7) = 2(f e ) + (U) - ^<R a • F a ) (2-19) 

i a 

3P N v = ^(R a • v a ur) = - Ytfi** • F «> ( 2 ' 2 °) 

a a 

= 2(f N ) - ^<R a • V a U) , (2-21) 

a 

since the external virials, YLii^i'^i^w) and J2a(^-a^ aU™ c ) , enable us to define the 
electron and nuclear pressures, respectively, on the base of the standard assumption: 
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Ylifci ■ VjC/ w ) = §Qy Pr ■ dS. 11 ^' 12 ^' 15 ^ For the volume parameter V to be a natural 
variable, the pressures, P e and Pn, must take a constant value at every point of the 
surface dV. In the above, Eq. (|2- 19|) is derived by the relation (|2 3|) . and the force 
on a- nucleus is given by F a = —V a U. The definition of the electron pressure ()2- 18j) 
can be used even when the nuclei behave as quantum particles without use of any 
approximation. 

At this stage, we can construct the following three systems of electron-nucleus 
mixtures by choosing the electron and nuclear wall potentials in different ways. 

(I) System-I, where the electrons and nuclei are confined in the volume V. In 
this system, the total pressure is given by 

3PV = 3(P e + PnW = 2{f c ) + 2(f N ) + (U) = 3PiV . (2-22) 

(II) System-II, where U^ e = and U™ c ^ 0. In this system, the electron 
pressure becomes zero because of (j2- 18p : 

3P e V c = = 2(f e ) + (U) - ^(R Q • F a ) . (2-23) 

a 

Therefore, the total pressure of this system can be attributed only to the nuclear 
part: 

3PV = 2(Tk) - ^<R Q • V a U) = 2(f N ) + ^(R Q • F a ) (2-24) 

a a 

= 2<f N ) + 2(f e ) + (U) = 3PnV , (2-25) 

At this point, it should be noted that the zero pressure relation (|2-23|) is used in the 
derivation of (f2 7 25|) from §TM$)- 

(III) System-Ill, where £/° uc = and U^ e = 0. This system represents nothing 
but an electron-nucleus mixture in the vacuum, where P e = and Pn = due to 
(|2T8p and (|2-20p . Consequently, the total pressure of this system becomes zero: 
P = P c + P N = 0. 

On the basis of the fundamental relations derived here in conjunction with 
System-I, II and III constructed in the above, we can make the following statements 
from (A) to (E): 

(A) : The electron pressure P e and the nuclear pressure Pn can be clearly defined for 
an electron- nucleus mixture by (|2- 19|) and (|2-2ip . respectively; the total pressure 
is a sum of the electron and nuclear pressures, as shown in (|2-22|) . 

(B) : The electron pressure is described by (|2-19|) . where the nuclear virial is sub- 
tracted. 

(C) : The pressure of System-I is the thermodynamic pressure of an electron-nucleus 
mixture. It has been shown that the pressure defined by the sum of (|2-19p and 
(|2-2ip becomes identical with the thermodynamic pressure determined by the 
volume derivative of free energy. 12 ) -17 ) 

(D) : The pressure Pi in System-I has an identical expression to Pn in System- 
II. Therefore, some investigators 1 )' 2 ^' 4 ) believe that System-II is identical to 
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System-I in the sense that both systems provide the same thermodynamic quan- 
tities (e.g., the pressure). However, the pressures, Pi and Pn, are equal with each 
other only when these systems are nearly equal to the system-Ill (P e = Pn = 0). 
It should be noticed that System-II is in general an inhomogeneous system, since 
this model involves a superatom in the extreme limit (V — > 0) with all nuclei 
compressed to become a point charge as this mathematical model. On the other 
hand, System-I is always homogeneous, since the wall effect can be replaced by 
periodic boundary conditions for this system. Therefore, System-I and System- 
II are different from each other in general, although (|2-22p and f)2-25|) have the 
identical expression; a physical quantity P e is not zero in the one and zero, in 
the other as Eq. (|2-23[) indicates. 
(E) : A metal in the vacuum is defined as an electron-nucleus mixture without the 
wall potentials to confine the electrons (U^ e = 0) and the nuclei (£7° uc = 0) in a 
finite volume. Therefore, by definition the thermodynamic pressure of a metal 
in the vacuum is zero, which is realized by the condition P c = and Pn = 0, 
as System-Ill exhibits: this results from the virial theorem for single particle, 
(|2T3p and (|2-14p . This fact makes a contrast to the usual treatment of a metal 
where the pressure of a metal becomes zero as a result of the cancellation in a 
sum: P = P e + P N = with P c < and P N > 0. 5 )' 6 > 



3. Electron and nuclear stress tensors 



In the present section, we derive formulae of stress tensors by the same arguments 
that were advanced in the preceding section treating the pressure in an electron- 
nucleus mixture. In the electron-nucleus mixture where the electrons and nuclei are 
confined in the same volume V (System-I), Nielsen and Martin [NM] 18 ' proved that 
the stress intrinsic to the system is given by a sum of the electron and nuclear parts: 



rpele ■ rpnuc 



(3-1) 



with the definitions of the electron and nuclear stress tensors by 



pele 
fiu 



E 



(3-2) 



- fj.u 



-E 



M 



■RtfVtJj) = -Nk B T6^-J2 (R^fA , (3-3) 



respectively. In the framework of the DF theory, the electron stress tensor is written 
in the form: 



pele 

- iiv 



V\xVv 



m 



o M r )p( r ))N , .„ drdr 



r — r 



'13 



+ J «;Mr + ^(P^P^)r 



(3-4) 
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where the bracket with N suffix indicates the configurational average of A({R Q }) 
defined by 

(A) N = ^ / dR N exp(-f3T c [n;{R a }])A({R a }) (3-5) 



with 



S = J dR N exp(-/3;F e [n; {R a } 



Here, the exchange-correlation stress tensor o~ xc is defined in the DF theory as 

in terms of the exchange-correlation part J- xc of the free-energy for electrons in the 
fixed nuclei, and /(e^) denotes the occupation probability for a state determined 
by the Kohn-Sham equation with its wave function (fi^. Also, the electrostatic energy 
£"es[ra] of the nucleus-electron mixture is written as 

^ . . e 2 f n(r)n(r') , , , A f Ze 2 n(r') , , 1 ^ Z 2 e 2 

E es \n] = — / , ; K ,/ drdr' - > / — y —jdr' + - > — (3-7) 

= E ee + E en + E nn = l [ Ppp^ldrdr' = E cs [n, n x ] (3-8) 
2 J |r — r'| 

with p(r) = — en(r) + Zeh\(r) and ni(r) = 5(r — R^). In the derivation of (|3-4p . 
we have used the relation: 

^ |p(r)p(rO (r ~^ ( ^~ r/) ^ r^ = -| n{v)r^ v ^dv -^R^E^ . (3-9) 

Therefore, by the combined use of (|3-3p and (|3-4j) . the total stress tensor is 
represented as 

-i | (p(r)p(rO) N (r ~ r0 ^ r^ + | (a^) N cfr , (3-10) 

which is essentially identical with the expression of the stress tensor given by NM if 
the nuclear kinetic term is neglected. 

On the other hand, in the system where only the nuclei are confined in the 
volume V (System II), the total stress tensor is given by 

V = T™ c = -Nk B T6^ - (R^Feu) N (3-11) 



which leads to the identical expression to (|3-10p . since the electron stress tensor T' 
becomes zero in this system. On the basis of System-II with the infinite nuclear 
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mass, Ziesche et al.^' and Dal Corso and Resta 19 ) have proposed to define the stress 
tensor in terms of the nuclear forces. In fact, Dal Corso and Resta have used the 
relation = in their derivation of (|3T0p from (|3-lip . In spite of (|3Tip leading 
to the same expression to (|3T0p . these formulae yield different results with each 
other for each system, since System-I is a homogeneous system while System-II is 
an inhomogeneous system in general. In a metallic system under high pressure, the 
electron stress tensor T™ 5 plays an important role compared with the nuclear stress 
tensor T™ c . 

§4. Definitions of electron pressure 

As mentioned in £j2j the virial theorem shows that the total pressure of an 
electron-nucleus mixture is given by a sum of the nuclear and electron pressures; the 
concept of the electron pressure is clearly defined by using the wall potential U^ e for 
the electrons in the form of (|2-19p , which is applicable to this mixture without use of 
the adiabatic approximation. On the other hand, many types of definitions for the 
electron pressure have been proposed up to date within the adiabatic approximation. 
At the first, the thermodynamic electron pressure is defined on the basis of the 
thermodynamic relation 10 ) 

, (4.1) 

TN e {n a } 



(I) P< 



' dV 



by introducing the electron free-energy J- e [n; {Ra}] for the electrons in the electron- 
nucleus mixture confined in the volume V: the nuclei are assumed as classical par- 
ticles, which produce an external potential for electrons in the system. This free 
energy, J- e [n; {R Q }], is defined by the electron part of the Hamiltonian for this mix- 
ture in the form 20 ) 

T c [n; {Ra}] = -fc B TlnTr c exp[-/3(f c + U)} , (4-2) 

where Tr e is the trace operator which is taken over a complete set of states of the 
electrons in the system; (|4-2p is the free energy of the electrons under the external 
potential caused by the nuclei fixed at {R a }. The DF theory provides an exact 
expression for the free energy of the electrons in this system in the form 

^ c DF [n; {R Q }] = F + E cs [n, + F xc [n, n{\ (4-3) 

with use of To, the free energy of a noninteracting electron gas, and the exchange- 
correlation contribution T xc to the free energy. 

In contrast to the above macroscopic definitions of the electron pressure, we can 
define the pressure in a microscopic way. In this spirit, Bader 21 ) -23 ) introduced the 
microscopic electron pressure tensor P e by treating an electron-gas in the fixed nuclei 
as an interacting many-body system in the form: 

(II) V • P c = F(r) , (4-4) 



where the force density F(r) is defined by use of the binary correlation function 
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r( 2 )(r,r') as 



F(r) = n(r)J^ V- 



Ze 2 



r-Rn 



r — r 



r, r 



(4-5) 



In a similar way, More 24 )' 25 ) defined the electron pressure tensor based on the fact 
that the electron gas is in hydrostatic equilibrium at each point in the electron- 
nucleus mixture with the pressure variations by electrical forces as is written in the 
following form: 



(III) 



VP e 



-ra(r)V 



N 

£| 

a=l 1 



Ze 2 



r-Rn 



+ 



e 2 n(r') 



r — r 



dr 



-n(r)V 



SE CS 
5n(r) 



(4-6) 



The same definition was introduced by Bartolotti and Parr 26 ) from a more general 
point of view. On the basis of the DF theory, the microscopic electron pressure tensor 
P e defined by (|4-6p is shown to be written explicitly by the sum of the electron kinetic 
tensor Pj^ F and the exchange-correlation pressure tensor Px C F : 27 ) 



P e ^ P° F + P X DF , (4-7) 
where P^ F is defined in terms of the wave functions </>j for the Kohn-Sham equation 



as 



pDF 



4m 



iVpVuti + ftVpVvfa - V^iV^* - V M #V„&] , (4-8) 



and the exchange-correlation contribution P^ F is defined in the DF theory as 



V • P° F = n(r)V 



5n(r) 



(4-9) 



In this definition, the thermodynamic pressure is given by the following surface 
integral over the surface of the system: 



3P P V = <b r ■ P e ■ dS . 
'dV 



which can be transformed in the volume-integral form: 27 ) 
3P e V = 2T S + E CS + [ trP^ dr - V R Q -F £ 



(4-10) 



(4-11) 



in terms of the kinetic-energy functional T s , where the operator tr means to take 
the sum of the diagonal terms of tensor. In the same scheme, another form of 
microscopic electron pressure tensor was introduced by Nielsen, Martin, Ziesche, 4 )' 18 ) 
and Folland 28 ) with use of the Maxwell pressure tensor P M 



(IV) 



~M 



PK F + Pxf + P M 



(4-12) 
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On the other hand, the electron pressure is determined by the formula based on 
the force theorem: 10 )' 29 )> 30 ) 



(V) 



+ 



3 DF 



dS, 



dv 



(4-13) 



where the volume derivative of the Kohn-Sham eigenvalue e, is performed under 
the self-consistent potential U se n being fixed. Prom a different point of view from 
the above pressure definitions, Janak 1 ) and Ziesche et alfi made a definition of the 
electron pressure with use of the total force on the nuclei in the form 



(VI) 



(4-14) 



In spite of different expressions, it can be shown that all definitions of the electron 
pressure, (I)-(V), except (VI), lead to the same pressure formula for the electrons 
confined in the volume V (see also Appendix [A]) as follows: 21 )' 23 )' 27 )' 10 ) 



3P e V = 2(f e ) e + (U)e ~ V(R Q ■ F a ) e = / r • P c 

~ JdV 



dS. 



(4-15) 



It should be noticed that the pressure defined by (|4-12)1 with use of the Maxwell 
pressure tensor can play the role of the electron pressure only under the condition 
that the relation, § QV r • P M • dS = 0, is fulfilled on the surface dV of the volume 
V. In fact, many investigators 31 ) -33 )' 30 ) use a surface-integral expression (|4-15p for 
calculating the electron pressure in solids. On the other hand, the definition (VI) 
leads to the relation 1 )' 4 ) 



3P,V = Y,R a F 



2(f c ) c + {U) c 



(4-16) 



with use of the adiabatic approximation and the Hellmann-Feynman theorem. Al- 
though the right-hand side of (|4-16p has an appropriate form of the virial theo- 
rem consisting of two times of the electron kinetic energy and the potential energy, 
Eq. (|4-16p cannot be taken as providing the electron pressure; this equation itself 
means P e = because of ([2-19P in Secl2j It is rather appropriate to consider that 
the pressure defined by (|4-14p provides the total pressure for another inhomogeneous 
system (System-II), that is, the electron-nucleus mixture where only the nuclei are 
confined in the volume V in the case where the nuclear kinetic pressure can be ne- 
glected (see also Appendix [Bj. In some works, 34 )' 35 ) the electron pressure is defined 
by the variation of the electron internal energy per unit cell with respect to the cell 
volume Q under the fixed entropy S and total electron number A^ e : 



(VII) 



Po 



dE. 



Q 



dQ 



(4-17) 



SN C 



This definition also leads to the relation 3P e V = 2(T e ) + (U), since the nuclear 
coordinates are not fixed in the ^-derivative of this definition (see Appendix|Aj) . Note 
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that the pressure (I4-17P is not the electron pressure, but provides the total pressure 
for an electron-nucleus mixture where the nuclear mass is taken to be infinite, as is 
shown in (|2 9|) . 

~M 

At this point, it is important to notice that the microscopic pressure tensor P 
defined by (|4- 12[) can provide a pressure formula for an arbitrary subspace with a 
volume fl in the system as follows: 



f : 

JdQ 



r • pf • dS 



/ trPf dr - V R Q -F a = <L r • P e • dS 

Jn a£Q JdQ 

2T s [n] n + E cs (f2) + / trP xc dr - V R a -F c 
Jo "To 



(4-18) 
(4-19) 



- M 



if P e is constant on the surface df2 of a volume O, and § dn r • P ■ dS = 0. Here, 
the electrostatic energy E es (f2) contained in the subspace f2 is defined as 



E CS (Q) = - [ p(r)r ■ V0(r)dr 
J a 



— I |E| 2 dr 



r • P M • dS 



(4-20) 
(4-21) 



on 



with use of the electric field E = — V(p and the electric potential 

,2 m /„A JV „ 2 



e0 r 



e n(r') , . Ze 
T\ dr 



r — r' 



a=l 



r-Rn 



5n(r) 



(4-22) 



In the above, it is important to recognize that the electrostatic energy E es (f2) con- 
tained in a subspace Q is not always equal to J Q |E| 2 dr/87r, as was defined in some 
works. 24 ) Also, the kinetic energy T s [n]fi involved in the volume fl is written from 
gSD as 25 )- 27 ) 



2T s [n] n = / trPg F dr 
J a 



£/(«) / 
2£/(«) / 

z J ft 



-/J 2 n2I 1 

V 2 4>i H 

2m 2m 



dr 



(4-23) 
(4-24) 



-fi 2 



ri 1 



V z d>idr+—I V 2 n(r)dr. (4-25) 
„ 2m Am J q K ' K ' 



In a similar way with use of the Maxwell pressure tensor, the electron part of 
the stress tensor is given for a subspace Q of the system by the formula 

T ele f2 = - [ Pfdr + VRfoF; = - I r opf. dS = - (f roP c .<iS, (4-26) 

Jn iea JdQ JdQ 



with the symbol o denoting a dyadic product of two vectors. 



14 



J. Chihara and M. Yamagiwa 



§5. Equation of states for liquid metals 



In ^21 we have proven that the pressure of an electron-nucleus mixture is ex- 
pressed by a sum of the electron and nuclear pressures, which are clearly defined 
in the virial forms. In the present section, we derive a virial pressure formula for 
simple liquid metals based on this fact, and the zero pressure condition is examined 
numerically for this pressure formula. 

5.1. Derivation of a virial pressure formula 

The standard virial pressure formula for a simple metal as an electron-ion mix- 
ture is written in a complicated form to calculate owing to the volume derivative 
of the effective ion-ion interaction in the formula. 36 ) In contrast to this standard 
expression, the virial theorem described in Sj2] provides the equation of states (EOS) 
in a simple form: 



3PV = 3(P C + P N )V = 2(T e ) - Y^{ii ■ ViU) + 2(f N ) - ^(R Q • V a U) (5-1) 

i a 

= 3P C V + 3Nk B T + ^ (R*F*)n > (5-2) 

a 

by regarding liquid metals as a mixture of electrons and nuclei. In a simple liquid 
metal, the interaction among ions with the Coulomb potential wfj(Q) can be ap- 
proximated as a sum of binary interaction vff(r), which is defined in the quantum 
hypernetted chain (QHNC) theory 37 ) by 

Pvff(Q) = (3vUQ) - — — i f * Q . (5-3) 

1 + Pn%v c ee {Q)[l - G{Q)\x u Q 

Here, C e \(Q) denotes the electron- ion direct correlation function, which plays the 
role of a pseudopotential w(Q) = —C c \{Q)/(3, and x%-> t ne density-density response 
function of a noninteracting electron gas, with G{Q) being the local-field correction. 
Therefore, with use of this two-body interaction, the EOS is written as follows: 

3PV = 3P e V + 3iV£;BT - N^n\ j g u (r)rVvff (r)dr . (5-4) 

At the next step, we proceed to derive an approximate representation of the 
electron pressure P e in (15-4)) on the basis of the spherical cell model. By taking a 
Wigner-Seitz sphere with a radius R, the electron pressure can be represented by 
the surface integral over the surface S in the form: 27 ) 

P « = ^ i r ' ( P K F + P ?c F ) ■ iS = P k (R) + Pxc(fl) - (5-5) 



3J? j s 

In the above, the kinetic pressure Pk(R) is given in units of hartree in the form: 33 ) 

P k(^) = ^I>(^-£)l<fe(#)l 



8tt 

ki 



R 2 



(5-6) 
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with 



1 



D ke = R— 



<t>ki dr 



(5-7) 



R 



and /x being the chemical potential. In the above, 4>ke( r ) denotes the radial part of 
the wave function for a state eke of angular momentum £ satisfying 



d 2 {r4> M ) 
dr 2 



+ 



2(e M - V eS ) 



£{£+!) 



(r 



0. 



(5- 



In the SJ model, 8 ) the whole crystal with dipole layer D at its surface is consid- 
ered as an assembly of Wigner-Seitz cells with dipole-layer surface. As a result, by 
taking the zero of a potential to be at the vacuum level, the effective potential V e s 
in the Wigner-Seitz sphere is written as 



V cS (r) =v cN + J v cc (r-r')n(r')dr' + Mxc (r) - D6(R - r) 



(5-9) 



with 9{x) being the step function and ^ X c( r ) = SJ-^ c /Sn(r). With a choice of the 
zero of the potential to be V c s(R~), the electron pressure at the surface S takes a 
form: 



P k (R) =-Dn(R) + ^L- e F )|fe(i?)| 5 



(D M -£)(D M +£+l) Q ■ 

^2 h le kl 



(5-10) 

In the above, ejL and e F indicate an energy level and Fermi energy with respect to 
this new origin, respectively, and the pressure at the surface is defined as Pk(R) = 

When the wave functions of electrons become plane waves and produce a con- 
stant density n near the boundary of the Wigner-Seitz sphere, the above equation 
leads to an approximate expression for the electron pressure as follows: 



Pe(n) 



-Dn + n 



■ dej 
dn ' 



(5-11) 



in terms of the energy ej per one electron in the jellium model. 



Table I. The nuclear pressures (GPa) for alkaline liquid metals in the vacuum calculated by the 
QHNC method and the Ashcroft model potential (Rb) . The nuclear kinetic and virial pressures 
are also shown together with the electron pressure given by PI" — n^dtj / dn\ n (the jellium 

model). The nuclear pressures calculated by Kumaravadivel 44 ' deviate significantly from the 
zero pressure condition compared with the QHNC results. 



' s 


rates T 


nuclear virial 


Pn 


JV 4) 


riodej/dn 




Li 


3.308 


0.2888 


-0.2871 


0.0017 


-1.7 


2.70 


Na 


4.046 


0.1253 


0.0362 


0.1615 


-0.24 


0.024 


Rb 


5.388 


0.0445 


0.0734 


0.1179 


-0.89 


-0.43 




5.388 


0.0445 


-0.0458 


-0.0013 




-0.43 
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In the jellium model, the electron pressure does not satisfy the condition P e = 
at equilibrium density no in the vacuum. In contrast to this model the condition 
P e = at equilibrium density uq is fulfilled in the stabilized jellium (SJ) model 7 )' 9 ) 
(ideal metal model 8 )' 9 )), where the finite jellium has an infinitesimally thin dipole 
layer on its surface. The relation (|5Tip provides the strength of the dipole layer at 
equilibrium [P e (no) = 0] of a metal in the vacuum: D(no) = [ndej /dn] no , which is 
identical with the result determined from the SJ model. In the S J model the strength 
of the dipole layer is assumed to be given by the relation: 



DM 



—D(no) . 
n 



(5-12) 



With use of this approximation the electron pressure in a liquid metal can be deter- 
mined finally by 



P e (n) 



n 



dej 


dej 




dn 


dn 


n _ 



(5-13) 



On the other hand, Nieminen and Hodge 38 ) showed that application of the re- 
lation (|5-6p to a free electron of plane waves leads to the expression for the pressure 
in terms of the internal chemical potential fl: 

P e {n) = n(/2 - ej) , (5-14) 

which can be also derived by the TF theory with the Weizsacher correction. 24 )' 26 )' 39 )' 40 ) 
At this point, Eq. (|5-14p can be rewritten in the form: 



n(fi - uj) + n 2 ^- (5-15) 
dn 

d(nej) /dn. By comparing this equation with (|5-lip we obtain the following 



with uj 
relation: 

D = fi 3 -fi, (5-16) 
which leads to the formula for work function u>: 41 ) w = D — uj = —p. with uj 
providing the internal chemical potential. From this equation we see that the relation 
—D/n = {fi — /Uj)/n = — dej/dn\ no is assumed in the SJ model. 

It has been already discussed about the problem that the electron pressure for a 
plasma does not become zero at equilibrium in the vacuum: in fact, in constructing 
the EOS for a plasma by the use of the ion-sphere model based on the TF method, 
a similar correction to (|5T3|) is introduced as a chemical bonding correction (Barnes 
correction) to make the electron pressure zero at equilibrium in the vacuum by More 
et aZ. 42 ) In the same spirit, we adopt (|5- 13|) as the electron pressure P e in ([5^ 
then, the total pressure for liquid metals is described in the final form: 



n 



dej 


dej 




dn 


dn 





+ nlk B T - i(n ) 2 J 9ll (r)vVvtf{r)dv . (5-17) 



At this point, it should be noticed that the following relation 

n\k B T = l -{nl) 2 J <m(r)rV<(r)dr 
must be satisfied at equilibrium in the vacuum, where Pn = 0. 



(5-18) 
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5.2. Numerical examination of the QHNC equation by the zero-pressure condition 

In this subsection, we exhibit how the zero-pressure condition, (15-18p . can be 
used as a check of numerical consistency in the calculation by applying it to some 
alkaline liquid metals at the normal pressure on the basis of the QHNC equation 43 ) 
as an example. Since we can regard liquid metals at the normal pressure as those in 
the vacuum, the pressures determined from (|5-17p must be zero in this circumstance: 
that is, Pn = while P c = by definition. In order to examine this relation we 
calculated the nuclear virial term and nuclear kinetic pressures from (|5-17|) . which 
are shown in Table U for Li, Na and Rb. 

In the case of Li, the nuclear virial term becomes -0.2871 GPa, which provides 
almost zero total pressure 0.0017 GPa by canceling out the nuclear kinetic pressure 
0.2888 GPa. Thus, the QHNC equation can be considered to give a reasonable result 
for Li to fulfill (|5T8p in comparison with -1.7 GPa of Kumaravadivel. 44 ) However, 
the QHNC equation does not provide even negative values of the nuclear virial term 
for the Na and Rb cases; as a result, the nuclear pressure remains to be a large 
positive value for each case. 

To see where this discrepancy comes from, we calculated the pressure of Rb 
from an effective ion-ion interaction determined by the Ashcroft model potential 
with the core parameter r c = 1.27A 36 ) on the basis of (|5-17p : the nuclear pressure 
becomes -0.0013 GPa, which satisfies the condition (|5-18|) to a considerable extent 
compared to other calculation 44 ) shown in Table HI Inversely, the condition (|5-18p 
can be used to determine the Ashcroft core parameter r c ; this condition generates 
r r = 1.275A for Rb case. In Fig. [H the effective ion-ion interactions for a Rb liquid 
metal obtained from the Ashcroft model with r c = 1.275 A are plotted together with 
that from the QHNC method. Although we find a visible difference between these 




i i i i i i i i i i i 

4 8 12 

r(A) 

Fig. 1. The effective ion-ion interactions of Rb (313K) determined by the QHNC equation and by 
the use of Ashcroft model potential (r c = 1.275), which are shown by the full and dotted curves, 
respectively. 
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two effective potentials in this scale, these potentials yield almost the same structure 
factors showing a good agreement with experiments as is displayed in Fig. [2j In 
contrast to the Ashcroft model potential, the pseudopotential — C e \{r)/ (5 determined 
from the QHNC equation has a complicate inner structure in the core region. As 
a consequence, the term J C^(r)dr [C^(r): the non-Coulomb part] which appears 
in (|5-18p is difficult to calculate with a sufficient accuracy; the discrepancy of the 
QHNC result for the zero-pressure condition is attributed to this problem in the Rb 
and Na cases in contrast to the case of a Li metal which has a simpler inner core 
structure in the pseudopotential C e \. It is well known that the inner-core structure 
of a pseudopotential is not significant to determine the effective ion-ion interaction, 
but is important to calculate the thermodynamical quantities, as the QHNC results 
indicate. It is worth noting that the condition (|5- 18j) provides a check of numerical 
consistency in the calculation for determining the effective ion-ion interaction, the 
structure factor and the thermodynamical quantities for a liquid metal in the vacuum. 

§6. Summary and discussion 

In general, the virial theorem (|2-10p teaches us that the partial pressure associ- 
ated with each component in a mixture can be uniquely defined by using the wall 
potential for each component, and the total pressure is represented as a sum of them: 
this definition is assured by the virial theorem for single particle, (|2-13p and (|2-14p . 
Also, the partial pressure can be actually measured as an osmotic pressure, for ex- 
ample. Therefore, in an electron-nucleus mixture, the electron and nuclear pressures 



o 

to 



0. 





1 | 1 1 1 1 | 1 1 




I Rb(313K) 




i .... i . . 







Q (A ) 



Fig. 2. The ion-ion structure factors Sn (Q) of Rb obtained by use of the effective ion-ion inter- 
actions based on the QHNC and Ashcroft models, which are displayed by the full and dotted 
curves, respectively. The full and white circles indicate the experimental results by the neu- 
tron 45 ' and X-ray 46 - 1 diffraction methods, respectively. The Ashcroft model with r c = 1.275A 
determined by the condition Pn = provides an indistinguishable structure factor from the 
QHNC result in a good agreement with the experiments. 
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are clearly and simply defined as the partial pressures owing to the virial theorem 
(|2-10p as is shown in £J21 and the thermodynamic pressure of this system is given 
by a sum of the electron and nuclear pressures. To be consistent with this fact, the 
electron pressure in this mixture is defined by 

3P C V = 2(f e ) - £<fi • V l U) = 2<T C ) N + (U C ) N - ^<R a • F a ) N , (6-1) 

i a 

which makes a contrast to another definition 1 )' 4 ) in terms of the forces on the nuclei 
3P C V = ^2(R a • F Q ) N = 2<T e ) N + (C4) N . (6-2) 

a 

It should be pointed out that the pressure represented by (|6-2p is neither the electron 
pressure nor the total pressure for this system, and that this pressure is associated 
with the total pressure of another inhomogeneous system (System-II as is named in 
52}, where only the nuclei are confined in the volume V and the electron can move 
over all space under the external potential caused by the confined nuclei. Similar 
remarks are made about the definition of the stress tensor for this mixture as is 
discussed in Sj3] 

In 21 we have enumerated several kinds of definitions for the electron pressure, 
which lead to the electron-pressure expression (|6-1|) . where the nuclear virial term is 
subtracted. It is important for the calculation of the total pressure by adding the 
nuclear pressure that we should recognize which definition is adopted to determine 
the electron pressure. For example 33 ) as was shown by (|4T5p . the electron pressure 
calculated by use of the surface integral expression, 3P e V = § d y r • P e • dS, is different 
from the pressure determined on the basis of the relation, 3P e V = 2T s [n]+E es + 
J trP^fdr, in that the latter contains a part of the nuclear pressure Yl a ^« ' ^ai 
which is approximated by the Madelung energy term: NCZ^ ut /4a as is given by 
Janak. 1 ) 

When the nuclei behave as classical particles, the partition function of an electron- 
nucleus mixture with the Hamiltonian (|2-ip is written exactly in the form 20 ) 



* 1 



N\h 3N 



dK N dP^exp 



(6-3) 



where the electron free energy J- e [n; {R a }] defined by (|4-2p plays a role of a many- 
body interaction among the nuclei: 10 ) 

F« = ~V Q T ?\tVN c = ^es| n ( r ) (6'4) 

Thus, the pressure of this system is obtained by performing the volume derivative 
of the free energy F = — ksT In without use of any approximation in the form 27 ) 



dF 

3PV = - 3V— 
oV 



3F tP) +3A^ B T+(^R a -F Q ) , (6-6) 



TN \ ®V I N 
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which is identical with (|5-2j) ; this relates the electron pressure determined in $4] under 
the adiabatic approximation to the total pressure and to that by the virial theorem 
(|5-ip owing to (|4-ip . That is, the relation, P = P c + P^, is proven without use of the 
wall potentials by (|6-6p in the adiabatic approximation. 10 ) When we approximate J- e 
involved in (|6-3p as 

F e [n; {Ra}} ~ Uq(V) + v n " R /3D • ( 6 - 7 ) 

and perform the volume derivative of F, we obtain the standard virial EOS 36 ) for 
simple liquid metals. On the other hand, the virial EOS (|5-17p is derived by applying 
the approximation (|6-7p only to the nuclear virial term of (|6-6p ; the first term of (|6-6p 
is treated appropriately as the electron pressure in contrast with the standard virial 
formula, 36 ^ in which derivation the approximation (|6-7p is applied to both terms 
involving T c in (16-6p . As a result, the electron pressure in the exact expression (|6-6p 
becomes to be represented by 

Po a + I(„.)3 / <m(r)?|£>A, (6-8) 

in the standard virial formula for a liquid metal. Therefore, the EOS (|5-17p is more 
accurate than the standard virial formula in spite of its simpler structure. 

A metal in the vacuum is defined as an electron-nucleus mixture without the 
wall potentials, which confine the electrons (U^ e ) and the nuclei (U™ c ) in a finite 
volume. From this definition, we can give another important remark that the zero 
pressure of an electron-nucleus mixture in the vacuum is realized by both electron 
and nuclear pressures being zero at the same time: this can also be verified directly 
from thermodynamics as described in Appendix [AJ In the adiabatic approximation, 
the electron pressure is given by thermodynamics in the form 

ZPeV ee - (3^ ^ c[n ^ R ° }] ) N = ({ 2T e + Ue -^K a .F a }) N 

= 2(f e ) - £(r ( • ViU) , (6-9) 

i 

which is zero in the vacuum due to (|213|) , since a metal in the vacuum is not confined 
by the wall. Similarly, = due to (|2-14p . Therefore, in the first-principles 
molecular dynamics (MD) performed on liquid metals in the vacuum, for example, 
the nuclear pressure should become zero (Pn = 0); this relation also can be used to 
check numerical consistency in the MD calculation. 

Appendix A 

Electron pressure determined by the DF theory 

The electron pressure can be expressed in terms of the internal energy of the 
interacting electrons in the presence of the external potential caused by the fixed 
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nuclei at {R a } as follows: 10 ) 



dT: 



DF 



dV 



8E V 



TN c {R a } 



dV 



SN c {R a } 



In the DF theory, this electron internal energy is written as 

E v [n; {R Q }] = T s v [n] + tf£[n; {R a }] + E^[n] . 
Here, the exchange-correlation part of the internal energy is defined by 



El c [n]=^M-T 



dF xc [n] 



dT 



Vn(r) 
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(A-l) 



(A-2) 



(A-3) 



and the kinetic-energy functional T^n] and the electrostatic energy E^ s [n; {R Q }] of 
this system are given by (|4-24j) and (|3-7p . respectively. With use of the scaled electron 
density distribution n\(r) = A 3 n(Ar), the electron pressure P e for this system can 
be determined by A derivative in the form: 



-3V 



dE v [n; {R Q }] 



dV 



A- 



dE x 



(A-4) 



A=l 



with 



E X = E v ^ [n x ; {R a }] = T// A3 [n x ] + E y J x3 [n x ; {R a }] + E^ [n x ] . (A-5) 

In the above, the scaled density distribution n\(r) keeps the total electron number 
as a constant in the A-variation: N c = J v n(r)dr = J V / X 3 n x (r)dr. With use of the 
following relations 

r/A S [n A ]=A 2 T s >], 
E^J x "[n x ; {R Q }] = XE^[n; {AR a }] 
and the exchange-correlation contribution 27 ) 



(A6) 
(A7) 



A 



dEXl X \n x ] 



A=l 



/ n(r)r • V 
Jv 



8n(r) 



dr + 



TV 



dV 



r-Pgf -dS 



we get finally the expression for the electron pressure in the form: 

dE x 



3P C V = A 



= 2TY[n] + El+ [ trP° F dr-^R a -F c 
1 Jv 

2T C + U C — ^2 R a • F a , 



trP° F dr, 
(A8) 

(A9) 
(A10) 



in terms of the true kinetic and potential energies given by T e = T s [n] + T xc and 
U e = E^ s + U xc , respectively. Here, the exchange-correlation contributions to the 
kinetic and potential energies, T xc and C/ xc , are defined by the relations 



trP^dr, 



E. 



v 

V 



(All) 
(A12) 
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On the other hand, in the adiabatic approximation, we have an average relation, 
(A) = ((y4.) e )N, where the bracket with a suffix e denotes the average over the electron 
states at a fixed nuclear configuration. Therefore, with use of this approximation 
and (|A-10p . the electron pressure f)2- 19|) defined by the virial theorem can be written 
as 

2(f e ) - £<fi • ViU) = ({2(T c ) c + (f7) c -^R a -F Q }) N (A-13) 

i a 

= ({ 2T e + U e - R« • Fa })n = 3<P c )nV , (A44) 

a 

which is identical with (|4-15p . It is important to note that the nuclei are fixed at the 
positions {R a } when the volume derivative is performed to determine the electron 
pressure. 

Appendix B 

Criticism on the Janak derivation of the pressure formula 



In the DF theory for the electron gas confined in the volume V [System- 1], the 
kinetic energy T s [n] is written as 

2T s [n} = - [ r-V-Pg F dr+/ r • Pg F • dS = [ trPjfdr (B-l) 
Jv JdV JV 

= [ n(r)r • W cff (r)eZr + / r • Pg F • dS , (B-2) 
JV JdV 

in terms of the kinetic pressure tensor Pg F defined by V • Pg F = -n(r)V£/"eff(r). 
Also, the exchange-correlation pressure tensor P xc , which is defined by the relation: 
V • P xc = n(r)V5E xc /5n(r), satisfies 

/ trP xc dr = - / r • V • P xc dr + I r • P xc • dS (B-3) 
Jv Jv JdV 

= - [ n(r)r • V^^-dr + I r • P xc • dS . (B-4) 
Jv ^ n ( r ) JdV 

Note that Eqs. (|B4]l and ([B~3l) are a general relation which must be satisfied for 
any tensor in the electron gas confined in the volume V. Nevertheless, the surface 
integrals in (|B-2p and (|B-4p are neglected in Janak's equations, 1 ^ (7) and (10) in the 
local-density approximation, respectively. 

Since the electrostatic energy E es , (|3-7p . fulfills the following relation: 10 ) 

f 5E 

E cs - VR Q -F Q = - / n(r)r- V— ^dr , (B-5) 
^— ' J v dn(r) 

with the force F a = —X7 a E es on the ath nucleus, and the effective external potential 
is defined by U e g (r) = 5[E es + E xc ]/5n(r), Janak has derived from (|B-2p and (|B-4p 

/ r • (Pg F + P xc ) • dS + V K a -F a = 2T s [n] + E cs + [ trP xc dr (B-6) 
JdV Jv 

= 2T C + U c , (B-7) 



Electron and nuclear pressures 
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with neglect of the surface-integral term in (IB-6P as mentioned before. If we put 
2>PV = ^R a -F Q in (|B-6P as Janak defined, we obtain 



3PV = 2T s [n] +E CS + [ trP xc dr - <f r • (Pg F + P xc ) • dS . 

JV JdV 



(B-8) 



Therefore, this pressure definition provides neither the electron pressure nor the total 
pressure of System-I in the Janak derivation. On the other hand, if we put in the 
above 

3P C V = I r • (Pg F + P xc ) • dS = 2T C + U c - V R Q -F Q = , (B-9) 
JdV 



the pressure given by (|B-8[) becomes the total pressure of System-II with infinite 
nuclear mass. In contrast with the above fact, Eq. (|B-6|) with the surface-integral 
term is identical with (|4- 15|) or (|2-9j) : the sum of the electron and nuclear pressures 
are equal to the total pressure of System-I with infinite nuclear mass. In reality, 
Janak's pressure formula for numerical calculation becomes that of System-I, since 
the surface-integral terms in (|B-2p and (|B-4|) , in comparison with his equation (7), 
are taken into account in his equation (23) at this stage. 



J. F. Janak, V. L. Moruzzi, and A. R. 
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